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Abstract 



We study the flux noise S$(oS) and finite frequency conductivity 01 (w) in 
two dimensional unfrustrated Josephson junction arrays (JJA's), by numer- 
ically solving the equations of the coupled overdamped resistively-shunted- 
junction model with Langevin noise. We find that S$(cv) oc w~ 3 / 2 at high fre- 
quencies ijj and flattens at low u, indicative of vortex diffusion, while o\ oc ui~ 2 
at sufficiently high oj. Both quantities show clear evidence of critical slowing 
down and possibly scaling behavior near the Kosterlitz-Thouless-Berezinskii 
(KTB) transition. The critical slowing down of 5$, but not its frequency 
dependence, is in agreement with recent experiments on Josephson junction 
arrays. 

PACS numbers: 74.50.+r, 74.40. +k, 64.60.Fr 



1 



Josephson junction arrays (JJA's) and thin-film superconductors are excellent model 
systems for studying vortex dynamics. At zero magnetic field, such systems are believed 
to undergo a Kosterlitz-Thouless-Berezinskii 1,2 (KTB) transition at a temperature Tktb- 
At temperatures below Tktb the vortices and anti- vortices are bound into pairs, whereas 
above T KTB these pairs start to unbind into unpaired vortices. Such a phase transition is 
expected to affect a variety of transport properties of such systems. 3 Measurements of both 
the IV characteristics and the inverse kinetic inductance consistent with the occurrence of a 
KTB transition have been reported in both thin superconducting films 4 and superconducting 
arrays. 5 

A particularly sensitive probe for such vortex dynamics is the study of flux noise. Con- 
ventional transport properties such as the IV characteristics, while sensitive to vortex dy- 
namics, are typically nonequilibrium measurements, requiring the application of an external 
current. By contrast, magnetic flux noise is typically measured at equilibrium, by placing a 
superconducting quantum interference device (SQUID) over a portion of the array or film. 
Such a measurement is sensitive to equilibrium fluctuations in the local vortex number of 
vortices within that area. 

A number of groups have studied flux noise in superconductors. Several measure- 
ments have been carried out in in high temperature superconducting films, including 
Bi 2 Sr 2 CaCu 2 8+ 5 6 and YBa 2 Cu 3 6 .95 J Recently, Shaw et al. s have done noise experi- 
ments on overdamped JJA's consisting of superconducting Nb islands in a Cu film, greatly 
extending some earlier measurements by Lerch et al. 9 These experiments yield a range of 
behavior for the spectral function S$(u) of the flux noise, that is, the frequency Fourier 
transform of the flux-flux correlation function. For example, YBa 2 Cu 3 6 .9 5 7 and JJA's 8 are 
found to have S$(u) oc uo^ 1 at "high" frequencies ("high," in this context, meaning greater 
than about 10-1000 Hz), while in Bi 2 Sr 2 CaCu 2 8 +(5 6 S$(u) oc cj~ 3 / 2 at similar frequencies. 

There have also been several theoretical studies of flux noise in such systems. Houlrik et 
al. 10 discussed the behavior of flux noise from a Coulomb gas analogy, and calculated S$(u) 
from a time-dependent Ginzburg-Landau (TDGL) model. At high frequencies, they found 

2 



S$(u) oc uo~ 2 . Gronbech- Jensen et al. 11 studied a JJA with a static magnetic field of 1/2 
flux quantum per plaquette, using the so-called resistively shunted junction (RSJ) model 
including self-capacitance. Their primary interest, however, was to find the voltage noise 
at finite external currents, with disorder in the islands' positions, rather than the flux noise 
itself. Recently, Wagenblast and Fazio 12 have studied the flux noise and scaling behavior 
using an XF-model with an assumed local damping for the phases. The local damping term 
corresponds to ohmic resistance shunts coupling each superconducting grain to the ground. 
The resulting flux noise was found to be white for low frequencies, varied as uj~ 2 at high 
frequencies, and as u^ 1 at intermediate frequencies. Very recently, Tiesinga et al. 13 have 
used both the coupled RSJ model and a TDGL model to study flux noise numerically over 
a relatively limited frequency range. They concluded that their TDGL results were closer 
to the experiment of Shaw et al. 8 than were the RSJ predictions. 

In this paper, we carry out extensive calculations of flux noise in an array of coupled 
overdamped Josephson junctions, using Langevin noise to simulate the effects of tempera- 
ture. Our model is similar to that of Tiesinga et al. 13 , but we study the real part of the 
frequency-dependent fluctuation conductivity <Ji(u) in addition to the vortex noise, and we 
calculate both over a considerably wider frequency regime. Our results for the flux noise 
show a clear signature of vortex diffusion above Tktb, i- e., S$(u) oc uj~ 3 / 2 above a cut-off 
frequency u> v (T) which approaches zero near T KTB . ai(u) is found also to have a character- 
istic frequency u) a {T) which approaches zero on either side of T KTB . 

The details of the RSJ model can be found in the literature. 14 The current through a 
junction between two superconducting islands i and j is assumed to consist of three contri- 
butions in parallel: a normal current Ir.^ = Vij/Rij through a resistance Rif, a Josephson 
current Is-ij = I c -,ij sm (%); an d a thermal noise current Here I c .^ is the critical current, 
Qij = &i — Qj is the phase difference across the junction, and = Vi — Vj is the voltage 
between islands % and j. The use of Kirchhoff's law for current conservation and of the 
Josephson relation = (h/2e)(d9ij/dt) leads to a set of coupled first-order nonlinear dif- 
ferential equations. We solve these coupled equations numerically by a standard algorithm 
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for square lattices of several sizes, as discussed, for example, by Chung et al. We assume 
no external current, and use periodic boundary conditions for an N x N square lattice of 
size N s = N x N with no disorder, i. e., I c -ij = I c and Rij = R for all and no external 
magnetic field. The time iteration is accomplished using a second-order Runge-Kutta pro- 
cedure with time intervals of order dt = 0.01r , where r = H/(2eRI c ) is the characteristic 
time of the problem. By means of the assumed Langevin dynamics, we can calculate not 
only time-dependent quantities, but also various equilibrium quantities 15 , computed as time 
averages. In the following, our results are presented in the "natural units" of the problem. 
Thus, the natural unit of time is To, and of frequency, ojq = l/r . Energy and temperature 
are given in units of KI c /(2e). 

Many of the physical observables of interest are spectral functions, that is, the Fourier 
transforms of various time correlation functions. For an observable 0(t), the spectral func- 
tion is defined as 



With this definition, So(^) has dimensions [O 2 ■ t], where t denotes time. 

Fig. 1 shows the time-averaged mean-square vortex density An z = (l/iV s ) X) r n z( r ) f° r 
zero magnetic field. Here n z (r) = or ±1 denotes the number of vortices in the r th plaquette, 
as conventionally defined. 16 An z becomes nonzero near T KTB , where the vortex-antivortex 
pairs begin to unbind. T KTB , as estimated from the vanishing of An z , is close to the accepted 
value of 0.9-0.95^/ c /(2e/c j B) for an infinite lattice. 17 

Another quantity of interest is the helicity modulus 7, 18 which measures stiffness against 
long-wavelength twists of 9, and is proportional to the superfluid density. Fig. 2 shows 7 and 
the integral of the fluctuation conductivity, 7 2 = (1/tt) ff^ ai(u)dw , where a\ = Recr. Both 
are calculated using standard expressions, 19 but from a time-average of the RSJ solutions. 
72 shows a characteristic peak near T KTB . The expected universal jump in 7 Tktb is 
somewhat broadened in our calculations, probably by finite-size effects. 

Fig. 3(a) shows o\{lo) itself, as calculated directly from the fluctuation-dissipation 
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FIG. 1. Time- averaged mean-square vortex density An z for the overdamped array at zero 
magnetic field, plotted as a function of temperature T for two different array sizes. 
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FIG. 2. Helicity modulus 7 and integrated fluctuation conductivity 72 for an overdamped array 
at zero magnetic field and two different array sizes, as calculated from time-averaged solutions to 
RSJ equations with Langevin noise. 
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FIG. 3. Fluctuation conductivity a\(uj) (eq. 2), plotted (a) versus frequency uj for several 
temperatures T, and (b) versus temperature at several uj. All quantities are plotted in natural 
units. The curves in (a) are vertically displaced; the dashed line has slope —2. 

theorem 20 for several different temperatures. In the classical limit (hw <C k B T) for this 
isotropic system, this theorem gives 

i f°° 

aM = ——Y, dtcos(ujt)(J x (r,t)J x (r',0)), (2) 

where J x = -f c sin(6>j — 9j) is the x-component of the supercurrent. For a fixed tempera- 
ture, <7i (a>) flattens at low frequency and falls off at high frequencies approximately as 1/u 2 . 
The slight upward convexity at high frequency here (and for the noise calculation below) is 
an artifact of the fast Fourier transform used to evaluate these quantities. Fig. 3(b) shows 
0i (uj, T) at several fixed cj's. At the lowest frequency (nominally uj — 0, but actually an 
average over several uj < 0.08), there is a strong peak at T ~ Tktb (somewhat masked by 
the log- log plot). For the other frequencies (all greater than u ), no peak is discernable near 
Tktb, indicating that the influence of the transition is suppressed at such high frequencies. 
By comparing these results with those of Fig. 2, we see that the peak in 7 2 is dominated by 
the low frequency regime of 0i(cj), i. e., uj < uj . 

Next, we turn to the flux noise S$(u). Rather than calculate this quantity, we instead 
calculate the somewhat simpler vortex number noise, S v (uj), which should behave similarly 
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FIG. 4. Same as Fig. 3, but for S v (lo). Dashed line has slope —1.5. 

in most cases. 8 S v (u) is the noise associated with N v = J2r n z( r ) 5 the vortex number enclosed 
in the area A spanned by the SQUID used to detect the flux noise. With periodic boundary 
conditions, N v would be identically zero if A were the entire array area. For a smaller 
A, N v fluctuates in time, giving rise to vortex noise. Our results are shown in Fig. 4 
for a square area A = I 2 equal to 1/4 of the array area. Fig. 4(a) shows our results 
for different temperatures as a function of frequency; (b), for different frequencies as a 
function of T. At low frequencies, S v becomes roughly frequency-independent, but at high 
frequencies S v (u) ~ c<j~ 3//2 , a dependence which is known to characterize diffusive behavior. 21 
The spectral function thus clearly shows the characteristic signature of vortex diffusion at 
temperatures above Tktb- 

Fig. 4(a) also suggests the same critical slowing down seen in the conductivity plots. 
That is, the S v (u) plots flatten out at a temperature-dependent crossover frequency u v (T) 
[inset of Fig. 5(a)] which vanishes as T — > Tktb from either side. Such behavior suggests 
that the effective vortex diffusion coefficient D — > as T — > T KTB . This follows from the 
approximate relation 21 D ps £ 2 u v which connects the D to the vortex noise spectrum in a 
fixed area I 2 . 

We now discuss more quantitatively the apparent critical slowing down seen in both our 
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FIG. 5. (a) S v (u,T) / S v (0,T) plotted against uj/u> v , where u v (T) is a crossover frequency ex- 
tracted from the data at each temperature and plotted in the inset, (b). a\{u>, T)/<7i(0, T) plotted 
against uj/lo^, where LO a (T) is plotted in the inset. 

flux noise and conductivity calculations. In the former case, Shaw et al. 8 have suggested a 
scaling form for S<s>(u) which may be written 



uS<s>(u) = F 



(3) 



where £ is the divergent correlation length characterizing the phase transition, and F is an 
appropriate scaling function (we ignore any dependence on other possible scaling variables 
such as £/£). The frequency is expected to vary as £~ 2 , where z(«2 experimentally 8 ) is 
a dynamical critical exponent. A simple scaling function which appears consistent with our 



calculations for S v (u) is 



F(x) = 



x 



(4) 



As a partial test of this form, we plot in Fig. 5(a) the ratio S v (cu)/S v (0) against u/u v , where 
uj v is a crossover frequency used as a fitting parameter to the numerical data. Clearly, all the 
plots fall atop one another, consistent with the scaling hypothesis, and are in good agreement 
with the form (4), which gives the diffusive form S v (u) oc u~ 3 / 2 at high frequencies. The 
crossover frequency cu v [inset of Fig. 5(a)] appears to go to zero near T KTB as required by 
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the scaling form, but our data are not adequate to test the expectation 8 uj v = oc £ 2 with 



£ oc exp byT B KT/(T — T B kt) • In short, our numerical results for S v at T > T KTB are 
consistent with diffusive vortex motion in the RSJ model, and a diffusion coefficient which 
vanishes continuously as T — > T KTB . 

The fluctuation conductivity <7i(u;) exhibits similar critical behavior, since it has the 
form of a Drude peak whose width goes continuously to zero as T — > T KTB . We may 
anticipate a scaling form <7i(u;) = o"o/(^ 2 + where a"o has critical behavior at Tktb] 
this form correctly gives the calculated 1/uo 2 high-frequency behavior. As a partial test of 
this hypothesis, we have plotted in Fig. 5(b), the ratio <Ti(u>, T)/<7i(0, T) versus uj/uj a for 
several temperatures, again using u a as a fitting parameter for each temperature [inset of 
Fig. 5(b)]. The calculated ai{u,T) all fall on the same universal curve both above and 
below T KTB (for clarity, we have plotted only the temperatures above T KTB ). Moreover, the 
crossover frequencies u a , like cu v , vanish smoothly as T — > Tktb from either side, though 
again our numerical data is inadequate to confirm the functional form u a oc with z ~ 2. 

Our results differ from the 1/uj behavior seen experimentally for S$(cj) in overdamped 
Josephson arrays 8 over several decades of frequency, although we do see indications of similar 
scaling behavior. Measurements on YBa 2 Cu30 6 . 9 5 7 show a similar 1/uj frequency regime, 
but as noted earlier, studies of very thin films of Bi 2 Sr 2 CaCu208+<5 6 also show a 
behavior in the same (100 Hz-lOKHz) frequency regime, which is also interpreted as evi- 
dence for vortex diffusion. Our results show clear evidence of such diffusive behavior in a 
numerical model for the dynamics of a Josephson array; the diffusing objects in our calcu- 
lations, as presumably in the Bi2Sr2CaCu20 8 +,5 films, are the thermally excited vortices and 
antivortices. The frequencies where we find diffusive behavior appear, however, to lie well 
above the diffusive regime in Bi 2 Sr 2 CaCu 2 8 +5 . 6 

It remains unclear why our calculations do not show a clear 1/uj regime, as seen in some 
experiments and reported in calculations using a local damping model. 13 One possibility 
is simply that we have not probed the noise spectrum to sufficiently low frequencies, or, 
possibly, in a sufficiently large area. In the experiments, the 1/uj regime occurs in the 
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KHz regime, and certainly well below MHz. By contrast, our calculations do not probe 
frequencies much below about 0.01u = 0.01 2eRI c /h. For typical array parameters, u may 
be in the range of MHz or higher, suggesting that it could be quite difficult to approach the 
1/uj regime in such numerical calculations. 

In conclusion, we have calculated both the vortex number noise S v (u) and the 
frequency-dependent conductivity ai(u) in an overdamped Josephson-junction array near 
the Kosterlitz-Thouless transition. The former exhibits a uj~ 3 ^ 2 behavior, characteristic of 
vortex diffusion, at high frequencies, and is flat at low frequencies. The latter has a Drude 
peak and a 1/uj 2 frequency dependence at high frequencies. Both quantities show clear ev- 
idence of critical slowing down (i. e., a vortex diffusion coefficient which goes to zero near 
Trtb) an d possible scaling behavior near the Kosterlitz-Thouless-Berezinskii transition, in 
agreement with experiment; but clear evidence of a 1/u regime for the vortex number noise 
is lacking. 

This work was supported by NSF Grant DMR94-02131 and DOE Grant DE-FG02-90 
ER45427 through the Midwest Superconductivity Consortium. Calculation were carried out 
using the SP2 at the Ohio Supercomputer Center. 
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